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Abstract 



X 

;-H I Phase transition in spherically symmetric collapse of a massless scalar field is 

studied in 4-d Einstein gravity. A class of exact solutions that show the evolution 
of a constant incoming energy fiux turned on at a point in the past null infinity are 
constructed to serve as an explicit example. The recently discovered phase transition 
[]I| in this system is manifest; above a threshold value of the incoming energy fiux, a 
black hole is dynamically formed and below that, the incoming fiux is refiected into 
the future null infinity. The critical exponent is evaluated and discussed using the 
solutions. 



Recent numerical investigations on gravitational collapse of a single massless 
scalar field in spherically symmetric Einstein gravity revealed a remarkable critical 
behavior near the onset of a black hole formation |jl]| [0. A generic one-parameter 
class of solutions S[p] decomposes into two phases depending on the magnitude of p 
that represents the strength of self-gravitational interaction of incoming scalar field. 
li p < p*, where p* is a threshold value, the gravitational collapse is followed by an 
explosion, reflecting back the incoming stress-energy flux into the future null infinity. 
For p > p*, a black hold is formed and its mass shows a universal scaling behavior 
Mbh — \p — p*\^ near the threshold p = p*. These approaches are based on a 
formulation of Christodoulou where this scattering process was studied as a Cauchy 
problem in general relativity |0]. In the context of 2-d dilaton gravity (CGHS model), 
a simplified theoretical model believed to capture many of the essential features of 4- 
d Einstein gravity, the explicit analytical understanding of the similar phenomenon 
is possible as the general (semi-)classical collapsing solutions are available and well 
understood 0]. To be specific, Strominger and Thorlacius found a similar universal 
critical behavior when Hawking radiation effect was taken into account [^]. The 
explicit analytical understanding of this critical behavior in 4-d Einstein gravity, 
however, has not yet been complete, partly due to the difficulty in solving 4-d 
Einstein gravity coupled with a single massless scalar field in a dynamical scattering 
situation. For example, the rigorous proof of the universality and the analytic 
calculation of the critical exponent remain to be studied. 

In this Letter, we make an attempt to analytically understand this critical be- 
havior. As a first step, we show that the constant incoming stress-energy fiux of 
scalar particles that was turned on at some point in the past null infinity can not be 
refiected forward into future null infinity if the magnitude of the incoming fiux ex- 
ceeds a certain threshold that plays a role oi p*. Rather, a black hole is dynamically 
formed. Below this threshold, the incoming fiux implodes through the origin and 
gets scattered forward into the future null infinity. This analysis supports the result 
of numerical studies, showing that this one-parameter class of solutions decomposes 
into the sub critical regime and the supercritical regime. The critical exponent that 
gives a scaling relation between the order paramer, a black hole mass in this case, 
and \p — p*\ can be calculated for these solutions in supercritical regime. Our ap- 
proach takes the classical back reaction of scalar matter fields on the space-time 
geometry into account exactly, since this critical phenomenon is the consequence of 
general relativity. Therefore, our result can not be observed in the leading order 
calculations of the weak field approximation. If we apply the weak field approxima- 



tion taking the Minkowskian metric as a fixed background, any unbounded amount 
of incoming flux will be imploded through the origin and all of it bounces forward 
into the future. 

First, we derive a class of exact dynamical solutions that can be utilized to 
construct an interesting physical situation. Assuming spherical symmetry, Einstein- 
Scalar action can be written as, 

/ = J d'xV^e-"^{R^'^ + 2g'^^d^<pdf,<P - 2e"f' - ^g-^djdfsf) (1) 

where we integrated out the angular coordinates using Gauss-Bonnet theorem. Here 

our convention is (H ) metric signature and the 4-d metric is given by ds"^ = 

Qapdx^'dx^ — e~'^'^{dO'^ + sin^ 9d(f'^). i?'-^^ and / represent the scalar curvature of the 
two dimensional longitudinal metric gaf3 and a massless scalar field, respectively. 
The gravitational constant G is IGnG = 1 in our unit. If we choose a conformal 
gauge for the longitudinal metric gap as gap = —e'^^~^'^dx^dx~, the action, Eq.(|l|), 
reduces to 

1 = J dx+dx-{And+d^p - ^ + ^d+fd-f) (2) 

along with two gauge constraints 

din - 2d±pd±n + ^n{d±ff = o 

where Q = e~'^'^. The general static solutions of the field equations obtained from 
the action (|^) under the gauge constraints can be readily found as follows; we can 
consistently reduce those partial differential equations into the coupled second order 
ordinary differential equations (ODE's) by assuming all functions depends on a single 
space- like coordinate x = x^x~ . The resulting ODE's can also be derived from an 
effective action 

dx{xtlp+\-^^-hlxh (3) 

and the original gauge constraints become 

VL-2ptl + -Qp = 0, 

where the dot represents taking a derivative with respect to x. The complete solu- 
tions to this effective action represent the general solutions of the static Einstein- 
Scalar fields under a particular choice of the conformal coordinates. To obtain these, 
we observe that the action (0) has three rigid continuous symmetries that allow us to 



construct three corresponding Noether charges, reducing the order of ODE's by one. 
The first symmetry is / — > / + a, which is clear as / field appears only through its 
derivative. The second symmetry is x -^ xe°' and p -^ p — a/2, which corresponds 
to the translation of the asymptotically flat spatial coordinate. The last symmetry 
transformation, x — > x^~^°', p —>■ p— ln(l+Q;)/4 — alnx/2, and Q — > Q{l + a), changes 
the action by a total derivative. This symmetry corresponds to the rescaling of the 
asymptotically fiat spatial coordinate. The Noether charges for these symmetries 
are constructed to be 

/o = xQf, (4) 

1-1 • 1 e^^ 

Co = x'^pn + -xn - -nx^f - -^x-^^, (5) 

and 

1 • 1 

xpQ — -xQ + -fl = cqIio-x + ci. (6) 

As these charges are conserved, /o, Cq and Ci are constants of integration independent 
of X. The gauge constraint, when combined with equations of motion from Eq.(|^), 
reduces to a condition cq = 0. We note that Eq.(|D can be readily solved for p and 
by putting this into Eq.(|^) we get a decoupled ODE. Once Q is solved in terms of 
X, p and / are easily determined as a function of x. If we take /o = 0, we recover 
static black hole solutions in Kruskal-Szeckers coordinates with ci = 4M^ where 
M is the mass of a black hole. For arbitrary value of /o, Ci > case solutions are 
identical to the solutions obtained by Janis et. al. using a technique to generate some 
Einstein-Scalar solutions from the vacuum solutions of the Einstein equations 0. 
Since ci = case contains no black hole as long as the static analysis is concerned 
and since the scalar field diverges logarithmically near the horizon in ci > cases, 
our calculation here is effectively a proof of no-scalar-hair theorem. The details of 
this consideration can be found in where complete static solutions of the more 
general action than Eq.([l|) are obtained. 

The Ci = cases are of particular importance as we can generalize it to dynamical 
situations we are interested in. The static solutions in this case are calculated to be 

p-2po 

n = ^(e^^o(ln(x/xo))^-4/o^), (7) 

p = —mil mx + pn, 

^ A 2 



„ , . Je^'^n + f^-fo ^ ^ ^ ^ ^ e^^°ln(x/xo)-2/o , , ^ 
/gSpof] + /2 + f^ e^po ln(a;/a;o) + 2/o 



where po, /i and xq are arbitrary constants. For simplicity, we take xq = 1, 
Po = and /i = for further discussions. Then, as i7 ^ oo, the behavior of / 
asymptotically approaches to — 2/o/r where the geometric radius r is defined to 
be V". In this limit, we find r -^ (Inx"*" + lnx~)/2 and the 4-d metric becomes 
ds"^ -^ —dx^dx~/{x^x~) — r^{d9'^ + sin^ Odip"^). After a conformal transformation 
Inx^ — *> x^, we find that the asymptotic space-time is a fiat Minkowskian. Addi- 
tionally, taking G ^ limit reproduces the same results. This limit is the same as 
the weak field approximation taking a Minkowskian metric as a fixed background 
metric. The wave equation for / in s-wave sector under this fixed background ge- 
ometry has general solutions of the form / = —2/o/r + constant. Thus, Eqs.(0) 
are the relativistic generalization of the point scalar charge solution. In a sense, 
Ci = limit is similar to taking an extremal limit of Reissner- Nordstrom black hole. 
This consideration suggests that the scalar charge /o, in some dynamic situations, 
can be a chiral field instead of being a strict constant, for /o is indeed an arbitrary 
chiral field in the weak field approximation. In general relativistic case, however, 
this simple generalization is not possible in general, since there can be a non-trivial 
corrections to Eqs.(0) of the order of d±fo. In case of asymptotically steady incom- 
ing and outgoing stress-energy flux, though, that simple generalization is possible. 
Forgetting about global boundary conditions, the result of this extension is 

n= -(Inx)^ - (A;+lnx+ - A;_ Inx" + go)^ (8) 

p = -lnfi lnx+ -ln(l + 4A;4.A;_), 

^ A 2 2 ^ ^ 

In X — 2(A;+ In x~*" — A;_ In x~ + go) N 

Inx + 2(/c+ lnx+ — /c_ lnx~ + go) 
where k± and go are constants. We can straightforwardly verify that these solutions 
satisfy field equations derived from Eq. (j|) and the corresponding gauge constraints. 
The additional constant term in the expression for p is the correction term originat- 
ing from d±f. This term turns out to be rather simple in this case where the charge 
/o has terms only up to linear terms in \nx^. The asymptotic stress-energy tensor 
averaged over the transversal sphere in a conformal coordinate system that becomes 
asymptotically flat near the past or future infinity is calculated to be T±± = Ak"^. 
Thus, 4fcj. are interpreted to be an incoming and an outgoing energy flux, respec- 
tively. The similar asymptotic analysis shows that go represents the background 
component of the scalar charge. Recently, go = and k^k_ = case of Eqs.(PD was 



reported in mathematics literature P] and named scale invariant solutions. For scat- 
tering situations in our consideration where incoming and outgoing flux can coexist, 
a slightly generalized version (§) proves to be useful. Additionally, the presence of go 
term enables us to consider the time evolution of the multiple square-type incoming 
energy pulses by successively gluing our solutions, unless a black hole is formed in 
an intermediate stage. To name a few other applications possible with our result, 
we can construct various scattering solutions, cosmological solutions and point par- 
ticle solutions with time-varying charge at the origin, depending on the boundary 
conditions and initial conditions. 

Using the results obtained so far, we construct the one-parameter class of exact 
solutions that reduce to / = {2klnx~^H{lnx^) — 2kln x~ H {— In x~)) /r as we take 
the leading order weak field approximation, i.e., taking G ^ limit. Here, H{x) 
denotes the usual step function. This approximate solution in a fixed Minkowskian 
background has a property that any unbounded amount of the incoming stress- 
energy fiux can be totally refiected off from the origin into the future null infinity. 
However, this picture qualitatively changes as we consider the exact solutions focus- 
ing on the space-time geometry change due to the stress-energy of the scalar field. 
The physical region of space-time in our consideration is specified by the requirement 
n > 0, since the angular coordinates should not have time-like signature. Thus, the 
natural boundary is fi = 0. The asymptotic incoming wave section of the weak 
field solutions is chosen as initial data and it represents the turn-on of the constant 
incoming energy fiux at a point in past null infinity. Under these conditions, the 
class of solutions parameterized by a constant A; > are found to be 

Q=-(u + vf I 
4 

fi = -(M + t,)2_A;V H (9) 

Q = j{u + vf - {kv - k^{k)uf III 

where we introduced f = Inx"*" and u = lnx~. Other functions can be read off from 
Eqs.(H). The specific form of k-{k) depends on the boundary condition at fi = 
in region HI. In our case, imposing a covariant version of the refiecting boundary 
condition yields k_{k) = (Vl + S/c — 16A;^ — l)/4. The region I, bounded by the 
past null infinity, v = and the origin u = —v, represents the Minkowski space 
before the turn-on of the constant incoming fiux. The region II, bounded by the 
future null infinity, u = and the past null infinity, represents the propagation of 
the incoming particles before any of them hits the boundary at the origin. Note 



that our coordinate choice sets the turn-on time of the incoming energy flux as 
V = 0. Finally, the region III, bounded by u = 0, the future null infinity and the 
origin u = —v{l — 2k) / {1 + 2k^{k)) , contains the scattering of the incoming particles 
off the origin and their further propagation toward the future null infinity. As the 
matter particles hit the origin, they linearly tilt it, adding some space-like component 
to its tangent line. li k < 1/2, the path of the origin remains time-like and the 
incoming particles are scattered forward into the future null infinity, just as in the 
subcritical regime of the numerical simulations. Especially, if A; <^ 1/2, the weak 
field approximation results are recovered and the geometry is almost Minkowskian. 
At k = 1/2, which can be interpreted as a phase transition point, i.e., the black 
hole formation threshold, the path of the origin, u = 0, becomes light-like for v > 0. 
li k > 1/2, the path of the origin becomes space-like and form a trapped region in 
space-time. In this case, the region III disappears and our solutions here become 
exactly the same as the one obtained in 0. As explained in detail in 0, the resulting 
space-time is a black hole with indefinitely increasing mass, for we do not turn off 
the constant incoming flux. Thus, this corresponds to the supercritical phase of this 
scattering system. As the numerical studies and the above considerations suggest, 
the order parameter of this system is the black hole mass. Then, the important 
physical quantity to compute, given our exact one-parameter class of solutions in 
supercritical regime, is the critical exponent. The geometric radius r = r^iv) of 
the apparent horizon of the dynamic black hole in supercritical case, determined by 
d,,r = 0, is calculated to be 



Jyl 



rH{v) = 2{k-^y/'{k+^-y/'kv. (10) 

The 1/(2G') = Svr times the value of this corresponds to the apparent mass Ma of 
the dynamic black hole. The linear dependence on v is understandable as it is the 
time duration between the turn-on of the incoming flux and the reference time v. 
We also note that the angular averaged incoming energy flux is Ak"^. Defining the 
transition point p* = 1/2 and p = k, we find that the critical exponent in this case 
is (5 = 1/2 in a scaling relation Ma — {p — P*Y ■ 

In the aforementioned numerical study [Q] , the incoming flux is a pulse type and 
because of this, the asymptotic out-region is a static black hole with asymptotically 
flat geometry. However, in our case, since we put inflnite mass into the black hole, 
the whole universe collapses to a future singularity, barring the existence of any 
realistic out-region. As a result, we should consider turning off the incoming flux 
at a flnite time f = fi > 0. Gluing a space-time with purely outgoing particles 



above v = vi is straightforward in subcritical case using Eqs.®. In supercritical 
case, however, the gluing is highly non-trivial, for we can not directly glue the 
static Schwarzschild solution at f = fi. For example, in 4-d Einstein gravity, gluing 
the Minkowski space to the Schwarzschild geometry requires an impulsive shock- 
wave type injection of matter particles 0. Thus, after we turn off the source, the 
geometry goes through a brief transient period to asymptotically settle down into a 
quasi-Schwarzschild geometry. During this process, the apparent horizon that was 
initially space-like at the turn-off time will settle down to a future null direction, 
slightly changing its geometric radius. The numerically obtained Mbh — {p — p*Y 
with 6 ~ 0.37 is very difficult to calculate, as the scaling relation between Mbh 
and Ma gets complicated through this process. The similar situation in the CGHS 
model is not as difficult as ours. The complicated transient behavior is absent in this 
case, due to the simplified dynamics of the model. Therefore, it would be possible 
to directly glue a static dilatonic black hole at f = f i and thereby getting a relation 
Ma ~ Mbh if we had considered CGHS model from the the outset. An interesting 
observation in this regard is that our 5 = 0.5 is the same as the critical exponent 
obtained by Strominger and Thorlacius in CGHS model 0. There is a reason for 



this connection as suggested in ||T0[. The pure gravity sector of CGHS model, other 
than being a target space effective action from string theory, can be considered as a 
leading order theory in the l/(i-expansion of the spherically symmetric rf- dimensional 
Einstein gravity. We can, therefore, adopt l/(i-expansion and consider the leading 
order behavior in the description of the complex transient process in 4-d Einstein 
gravity. As a zeroth order approximation, we glue a CGHS black hole directly to our 
solutions to deduce the approximate scaling relation Ma — Mbh- Thus, the leading 
order approximation of the exact critical exponent for finite pulse-type incoming 
energy fiux is now calculated to be 0.5. Since the next order correction to the critical 
exponent is expected to be an order of 1/d = 0.25 and the numerically calculated 
value is about 0.37, our leading order value, 0.5, seems plausible. Generalizing this 
consideration, it is conceivable that the exact critical exponent for finite pulse-type 
incoming energy fiux in rf- dimensional spherically symmetric Einstein gravity will 
interpolate 5(4) ^ 0.37 in 4-d Einstein gravity and 6{oo) = 0.5 in CGHS model as 
a continuous function 6{d) of the space-time dimensionality d. In many other cases 
of phase transitions, the critical exponent rapidly saturates into a limiting value at 
d = oo. It will be an interesting exercise to verify this conjecture and, additionally, 
develop a systematic perturbation theory with a dimensionless expansion parameter 
1/d to tackle other difficult problems in 4-dimensional gravity. 
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